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ABSTRACT 

Excited  state  proton  transfer  dynamics  are  reported  for  the  l-naphthol(NH2)n 
cluster  system  for  n  =  3  and  4.  Picosecond  time  and  mass  resolved  pump  (Sj  <-  Sq)  probe 
(I  <—  Sj)  experiments  demonstrate  the  following  results:  1.  excited  state  proton  transfer 

occurs  for  n  =  3  and  4  clusters  only;  2.  for  n  =  5  clusters  the  proton  is  transferred  in  the 

ground  state  and  for  n  =  2  clusters  no  proton  transfer  can  be  observed;  3.  the  proton 
transfer  time  in  the  n  =  3  cluster  at  the  Oq  position  is  ca.  60  ps;  4.  this  time  is  reduced  to  ca. 

40  ps  and  ca.  10  ps  for  800  cm'^  and  1400  cm‘*  of  vibrational  energy  in  Sj,  respectively; 

5.  for  the  n  =  4  clusters  these  times  become  approximately  70,  70  and  30  ps,  for  0,  800 
and  1400  cm'*  of  vibrational  energy  in  Sj,  respectively;  6.  both  n  =  3  and  4  clusters  exhibit 

a  second  low  amplitude  decay  component,  which  is  about  an  order  of  magnimde  slower 
than  the  initial  decay;  and  7.  l-naphthol-djCNDj)^^  clusters  have  a  greatly  reduced  rate 

constant  for  the  excited  state  proton  transfer  dynamics.  These  observations  are  well  fit  and 
explained  by  a  simple  statistical/barrier  penetration  model  involving  proton  tunneling  and 
the  effect  of  van  der  Waals  vibrations  on  the  height  and  width  of  the  barrier  to  proton 
transfer. 


I. 


INTRODUCTION 


Proton  transfer  is  a  central  mechanistic  and  kinetic  step  in  a  large  number  of 
important  reactions.  Reactions  involving  proton  transfer  are  strongly  affected  by  the 
surrounding  solvent:  several  studies  of  proton  transfer  mechanisms  and  dynamics  in  bulk 
condensed  phases  have  appeared  encompassing  solvent,  pH,  and  isotope  effects.*’^  In 
these  studies  the  effects  observed  can  be  related  to  the  bulk  properties  of  the  solvent.  The 
individual  effects  of  solvent  molecules  and  solute-solvent  geometry  can  also  be  explored  in 
the  condensed  phases  through  the  study  of  matrix  isolated  clusters.^  Gas  phase  studies  of 
clusters  and  cluster  chemistry  are  also  now  appearing,^'^  which  can  provide  additional 
microscopic  details  of  the  solute-solvent  structure  and  the  affect  of  solvent  on  solute- 
solvent  chemistry  (e.g.,  proton  transfer).  These  latter  studies  are  of  both  a  static^'^  and 
dynamic  ’  nature.  Two  very  sensitive  tests  of  the  mechanism  of  a  proton  transfer  reaction 
are  the  hydrogen/deuterium  isotope  effect  on  the  reaction  rate  constant  and  the  effect  of 
vibrational  excitation  within  the  solute-solvent  cluster  on  the  reaction  rate  constant.  Proton 
transfer  reactions  in  gas  phase  clusters  are  not  well  explmed  partly  because  these  latter  two 
experiments  have  not  been  performed  and  partly  due  to  other  experimental  dif^culties. 
Studies  to  date  suggest  that  proton  transfer  for  phenol  and  1-naphthol  clustered  \»ith  water 
and  ammonia  solvents  is  a  function  of  cluster  size  and  composition.^'^ 

This  report  is  concerned  with  proton  transfer  in  1-naphthol/ammonia  clusters;  in 
particular,  the  effects  of  cluster  size,  vibrational  state,  and  hydrogen/deuterium  substitution 
(i.e.,  l-naphthol-djfNDj)^  vs.  l-naphthol-hjfNHj)^^)  on  the  excited  state  proton  transfer 

dynamics  are  observed  and  modeled.  Excited  state  proton  transfer  is  found  to  take  place  in 
these  systems  fw  clusters  with  three  and  four  ammonia  solvent  molecules.  The  transfer  rate 
increases  for  increased  Sj  vibrational  excitation,  decreases  for  the  deuterium  substituted 

clusters,  and  is  independent  of  the  wavelength  of  probe  laser  used  to  ionize  the  cluster  for 
detection.  A  simple  model  involving  proton  tunneling  is  proposed  which  explains  these 
observations  nearly  quantitatively. 


II.  EXPERIMENTAL  PROCEDURES 


Naphthol/ammonia  clusters  are  formed  in  a  seeded  helium  supersonic  jet  expansion. 
1-Naphthol  is  purchased  from  Kodak  and  used  without  purification.  For  isotope  effect 
studies  the  hydroxyl  proton  of  1-naphthol  is  exchanged  with  deuterium  by  dissolving 
naphthol  in  anhydrous  ether  and  washing  the  solution  several  times  with  D2O.  The  ether 

and  water  layers  arc  separated  from  one  another  and  the  ether  is  evapora^  under  vacuum. 
1-Naphthol-OD  (l-naphthol-d^)  crystals  are  recovered  from  this  process. 

1 -Naphthol  is  placed  in  the  sample  compartment  of  a  pulsed  nozzle  and  heated  to 
roughly  80®C.  The  helium  expansion  gas  is  maintained  at  55  psig  and  contains  0.05  to  0.1% 
ammonia.  The  molecular  beam  enters  the  ionization  region  of  a  resonance  enhanced,  two- 
photon  ionization,  mass  spectrometer.^®  The  clusters  formed  in  the  expansion  are  detected 
by  a  two-color  pump/probe  ionization  scheme  involving  the  transitions  I  <-  Sj  <-  Sq.  Time 

resolved  signals  for  the  Sj  state  proton  transfer  reaction  are  obtained  by  varying  the  time 

delay  between  the  pump  and  probe  beams. 

The  laser  systems  used  for  this  study  is  shown  in  Figure  1.  A  mode-locked, 
frequency-doubled,  cw,  Nd:YAG  laser  (Spectra  Physics  model  38(X)S)  synchronously 
pumps  two  dye  lasers  (Spectra  Physics  model  35(X)).  Approximately  10%  of  the 
remaining  1064  nm  beam  is  used  to  seed  a  NdrYAG  regenerative  amplifier  (Spectra 
Physics  model  GCR-3A),  the  output  of  which  when  frequency  doubled  provides  up  to  40 
mJ  of  532  nm  light  in  a  50  ps  pulse  at  a  repetition  frequency  of  10  Hz.  This  pulse  is  split 
50/50  and  pumps  two  home-built  4-stage  dye  amplifier  trains.  The  outputs  of  the  amplifier 
trains  are  tunable  pulses  of  ~5  to  8  ps  duration  with  energies  of  1-2  mJ.  These  pulses  are 
frequency  doubled  in  type  I  KDP  crystals  to  give  tunable  UV  radiation  with  an  intensity  of 
a  few  hundred  microjoules/pulse.  The  UV  pulses  can  be  Stokes  Raman  shifted  with  nearly 
30%  efficiency  in  CH^  to  extend  the  accessible  wavelength  range  of  the  system.  Either  one 


of  the  beams  can  be  delayed  with  respect  to  the  other  using  a  computer  controlled,  stepper 
motor  driven  linear  motion  stage. 

Data  acquisition  is  accomplished  by  setting  an  integrator  gate  on  the  appropriate 
time  of  flight  for  Jie  ion  of  interest.  The  output  of  the  integrator  is  sent  to  an  analog-to- 
digital  converter  interfaced  to  the  computer  which  ctmtrols  the  delay  stage. 

The  excitation  laser  uses  kiton  red  dye  to  give  excitation  energies,  after  frequency 
doubling,  in  the  region  from  30,800  cm'^  to  32,500  cm'^  The  amplifier  train  employs 
either  DCM  or  R640  as  the  gain  medium.  The  ionization  laser  contains  DCM  and  is 
amplified  by  either  DCM  or  LDS  698  and  frequency  doubled.  The  ionization  energy  used 
for  most  experiments  is  in  the  range  from  28,500  to  29,500  cm'^  The  ionization  laser 

wavelength  can  be  Stokes  shifted  (2914  cm**)  by  focussing  the  beam  intt)  a  cell  containing 
120  psi  of  CH^,  permitting  formation  of  the  cluster  ion  with  lower  vibrational  energy.  The 

excitation  laser  pulse  energy  is  kept  below  20  microjoules^ulse  for  all  experiments  while 
the  ionization  laser  pulse  energy  is  maintained  at  100-200  microjoules/pulse.  The  pulse 
width  of  the  amplified  dye  beams  can  be  determined  by  autocorrelating  or  crosscorrelating 
the  amplified  beams  using  the  delay  stage  aiKi  a  KDP  doubling  crystal.  The  dye  amplifiers 
are  run  at  high  gain  for  experiments  requiring  Raman  shifted  beams  generating  350-400 
microjoules/pulse  of  frequency  doubled  dye  laser  energy  resulting  in  90-100 
microjoules/pulse  in  the  Raman  shifted  beam.  Such  high  gain  in  the  amplifier  chain 
broadens  the  ioruzation  laser  pulse  by  a  factor  of  ~  3  to  ca.  20  ps. 

III.  RESULTS 

A.  Spectroscopy 

The  steady  state  excitation  and  ionization  spectra  of  1-naphthol/ammonia 
clusters  have  been  reported.^’^*  These  spectra  provide  the  necessary  excitation  and 
ionization  energies  for  the  clusters  of  interest  as  a  function  of  cluster  size.  The 
l-naphthoKNHj)!^  clusters  display  sharp  transitions  that  are  red  shifted  from  the  bare 


molecule  1-naphthol  origin  (  31,458  cm'^)  by  240  and  350  cm'^  respectively.  The  n  =  3 

clusters  possess  both  sharp  and  broad  features  near  their  origin  transitions  both  of  which 

are  further  red  shifted  from  the  bare  molecule  origin.  The  n  =  4  clusters  have  only  broad 

spectra  still  further  red  shifted  from  the  1-naphthol  origin.  The  n  =  5  cluster  shows  an 
anomalously  large  (1000  cm'^)  red  shift  relative  the  bare  molecule  (^. 

The  ion  appearance  potentials  for  the  n  =  1  and  2  ammonia  cluster  ions  are 
red  shifted  ~  5(XX)  cm*^  from  that  of  the  bare  1-naphthol  molecule  (62,600  cm'*).  The  ion 
appearance  potentials  for  the  different  geometries  of  the  n  =  3  clusters  are  red  shifted  from 
that  of  1-naphthol  by  ~  5000  cm'*  for  the  cluster  that  does  not  undergo  excited  states 
proton  transfer  and  -  6600  cm'*  for  the  cluster  that  does  undergo  proton  transfer.^^ 
Clusters  with  n  >  4  all  have  the  same  ionization  appearance  potential  of  ~  55,900  cm'*. 
The  threshold  for  proton  transfer  as  a  function  of  cluster  size  for  the  1-naphthol/ammonia 
system  has  thus  been  suggested  to  be  that  geometry  of  1-naphthol (NHj)^  with  the  lower 

(like  n  =  4,  ...)  ionization  appearance  potential.  The  previous  time  resolved  studies^ 
further  substantiate  this  assignment 

1-Naphthol  is  known  to  have  two  rotomers;  syn  and  anti  orientations  of  the 
hydroxyl  group.**  The  anti  conformer  is  the  mote  stable  and  is  the  only  one  which  has 

significant  population  in  our  supersonic  jet  expansion.^ 

B .  Dynamics  -  l-naphtholfNIlj)^,  n  =  2-5  Origin 

Time  and  mass  resolved  excitation  (1  +  1,I<— Sj<—  Sq)  data  are  collected 
on  the  n  =  2-5  clusters  at  the  energies  of  their  respective  Sj  «—  Sq  origin  transitions.  Only 

the  n  =  3  and  4  clusters  exhibit  any  time  dependent  results.  The  time  resolved  results  are 
independent  of  ionization  laser  wavelength  within  the  range  from  near  the  ion  appearance 
threshold  (25,(X)0  plus  500  cm'*)  to  5000  cm'*  above  the  threshold. 

Figure  2  shows  a  typical  1  +  t  time  and  mass  resolved  ionization  signal  for 
l-naphtholfNHj)^.  The  data  are  fit  to  a  bi-exponential  curve  with  tj  =  60  and  X2  =  500  ps. 


and  relative  amplitudes  of  =  0.3.  The  instrument  response  function  (~  16  ps  fwhm) 

is  obtained  from  the  rise  of  the  n  =  1  and  2  cluster  signals.  The  decay  time  constants  are  fit 
with  a  Levenberg-Margnardt  algorithm.  The  observed  time  dependence  is  assigned  to  the 
proton  transfer  process  based  on  its  cluster  size  dependence  and,  as  discussed  below, 
isotope  effect.  The  possibility  of  cluster  fragmentation  has  been  eliminated  by  careful 
searches  for  fragmentation  products.  The  observed  decay  must  therefore  be  due  to  a 

change  in  ionization  cross  section  as  the  proton  is  transferred  from  1-naphthol  to  ammonia 
in  the  Sj  state.  The  1-naphthol  (NHj)^  cluster  decay  times  at  the  origin  are  similar: 

Tj  ~  70  ps  and  ~  800  ps,  with  a  relative  intensity  ratio  of  A^Aj  ~  1.0.  These  data  are 
summarized  in  Table  I.  As  in  the  n  =  3  case,  these  decay  time  constants  arc  independent  of 

ionization  energy  from  500  to  4000  cm"^  above  the  ionization  appearance  threshold. 

The  analysis  of  the  proton  transfer  dynamics  of  l-naphthoKammonia)^  ^ 

discussed  below  will  focus  on  the  proton  transfer  decay  time  t  j.  The  second  long  decay 
time  Xj  can  be  attributed  to  either  a  cluster  reorganization  (which  changes  the  I «-  S  j  cross 

section)  subsequent  to  the  proton  transfer  event  or  a  proton  transfer  associated  with  a 
cluster  of  different  geometry. 

The  proton  transfer  time  (Xj  ~  60  ps)  for  the  gas  phase  l-naphthoUNIlj)^  cluster  is 

a  factor  of  2  to  3  slower  than  the  time  observed  for  the  argon  matrix  isolated  cluster  (20  - 
25  ps).^*’  The  faster  proton  transfer  rate  in  the  matrix  system  can  probably  be  attributed  to 
polarizability  of  the  surrounding  environment.  Matrix  polarizability  can  have  two  effects 
on  the  proton  transfer  reaction  rate;  thermodynamic,  in  which  the  matrix  stabilizes  the  ion 
pair  products;  and  kinetic,  in  which  the  matrix  reduces  the  height  of  the  barrier  to  proton 
transfer. 

C.  Dynamics  -  l-naphthol(NH3)jj,  n  =  2-5  Vibronic  Transitions 

Excitation  of  the  1-naphthol  molecular  vibrations  in  the  Sj  state  of  the  van 

der  Waals  cluster  is  observed  to  increase  the  proton  transfer  rate  even  though  the  0-H 
stretch  motion  itself  is  not  excited,  and  is  not  energetically  accessible.  Specifically,  if 


~  900  cm  of  vibrational  energy  is  placed  in  the  n  =  3  cluster  Sj  state  the  signal  decay 
times  become  tj  ~  40  ps  and  X2  -  500  ps.  At  ~  1400  cm'^  of  vibrational  energy  for  the 
n  =  3  cluster  in  Sj,  the  decay  times  drop  to  Xj  ~  12  ps  and  X2  ~  120  ps.  For  the  n  =  4 
cluster  the  decay  times  at  ~  1250  cm'^  of  Sj  vibrational  energy  are  Xj  ~  70  ps  and  X2  ~  600 
ps  and  at  ~  1700  cm'^  of  S  j  vibrational  energy  the  decay  times  become  Xj  ~  30  ps  and  X2 

~  1000  ps.  These  data  are  also  summarized  in  Table  I.  The  n  =  4  cluster  signals  are  much 

weaker  and  have  much  larger  error  limits  than  the  n  =  3  signals.  Nonetheless,  they  still 
show  that  the  n  =  4  cluster  has  a  reduced  X|  time  with  increasing  vibrational  energy  and  that 

the  overall  effect  of  vibrational  energy  is  less  than  that  observed  in  the  n  =  3  cluster. 

Even  with  up  to  1700  cm**  of  vibrational  energy  in  Sj  and  up  to  4000  cm'* 
of  vibrational  energy  in  the  cluster  ions  no  fragmentation  of  l-naphthol(NH2)3  ^  is 
observed:  we  do  not  find  (NH3)|jjH‘‘'  signals  nor  is  any  time  dependence  found  in 
l-naphtho^NHjlj  2  mass  channels. 

D.  Dynamics  -  l-naphthol-djfNDj)^ 

Perhaps  the  most  convincing  demonstration  that  the  observed  time 
dependent  cluster  ionization  signals  for  l-naphtholCNH^)^  ^  are  indeed  due  to  excited 

state  proton  transfer  is  a  substantial  deuterium  isotope  effect  on  the  observed  rates. 
Moreover,  if  the  changes  in  decay  times  are  large,  one  can  expect  tunneling  to  be  a  major 
contribution  to  the  proton  transfer  mechanism. 

The  observed  decay  times  Xj  and  Xj  for  the  deuterated  n  =  3  cluster  at  the 
origin  and  at  1400  cm**  of  vibrational  energy  in  Sj  are  given  in  Table  H.  For  the 
l-naphthol-djfND^)^  cluster  excited  to  the  Sj  <—  Sq  Oq,  virtually  no  proton  transfer  is 

observed.  At  1400  cm**  of  vibrational  energy  the  decay  time  Xj  of  the  signal  decreases  to 
75  ps,  roughly  a  factor  of  six  larger  than  observed  in  the  l-naphtholCNHj)^  cluster. 

These  results  corroborate  the  proton  transfer  interpretation  of  the  observed 
Sj  l-naphtholCNHj)^  ^  decay  times  and  additionally  suggest  a  mechanism  for  this  reaction 


which  involves  proton  tunneling  through  a  barrier.  Similar  results  and  conclusions  obtain 
0 

for  the  Oq  studies  of  matrix  isolated  l-naphthoKammonia)^  clusters. 
rV.  DISCUSSION 

In  this  section  we  present  a  model  for  calculating  the  proton  transfer  rate  constants 

for  van  der  Waals  clusters.  The  model  employed  consists  of  three  essential  features:  1.  the 

untransferred  and  the  transferred  structures  are  separated  from  one  another  by  a  potential 

energy  barrier  which  can  be  characterized  by  a  width  and  a  height;^’  2.  the  barrier 

width  and  height  are  modulated  by  vibrational  excitation  of  the  intermolecular  cluster 
IS  1 

modes;  and  3.  vibrational  energy  is  distributed  statistically  among  the  vibrational  (van 
der  Waals)  modes. 

Two  vibrational  coordinates  take  on  special  signiflcance  in  this  model.  First,  the 
reaction  coordinate  is  taken  to  be  the  OH  stretch.  The  hydroxyl  group  is  hydrogen  bonded 
to  the  ammonia  nitrogen,  and  the  barrier  to  proton  transfer  is  taken  to  be  along  this 
O-H  -N  reaction  coordinate.  Second,  the  O— N  distance  in  this  hydrogen  bond  is  of 
central  importance  in  determining  the  height  and  width  of  the  proton  transfer  barrier.  This 
distance  can  be  modulated  by  excitation  of  the  "diatomic"  stretch  mode,  for  which  the  1- 
naphthol  and  (NH^)^  moieties  are  taken  to  be  the  "atoms".  Excitation  of  the  O  -N  stretch 

can  dramatically  affect  the  proton  tunneling  rate.^^'^*  Numerical  application  of  this  model 
requires  calculation  of  the  proton  tunneling  rate  for  each  quantum  state  of  the  diatomic 
stretch  mode  and  calculation  of  the  probability  of  occupation  of  each  quantum  state.  These 
calculations  are  discussed  below.  Tunneling  rates  can  be  calculated  as  a  function  of  the 
O  -N  separation  based  upon  the  WKB  approximation  for  particle  penetration  through  a 
barrier  of  assumed  functional  form. 

The  model  potential  surface  for  this  system  is  presented  in  Figure  3  ir  which  the 
untransferred  proton  is  bound  in  a  well  and  the  transferred  proton  is  represented  as  a  free 
particle.  The  surface  in  Figure  3  thus  represents  the  hydrogen  stretch  coordinate.  The 


proton  transfer  reaction  is  exothermic  and  the  products  are  expected  to  have  a  high  density 
of  states:  this  justifies  both  the  "free  particle"  product  description  and  the  barrier 
penetration  limit  Under  these  conditions,  the  tunneling  probability  or  rate  constant  for  the 

proton  transfer  reaction  reduces  to  the  well  known  expression*"* 

k  =  vexp[-|-(2mE^)*/2]  (1) 

in  which  v  is  the  zero  point  frequency  for  the  proton  in  the  potential  well,  is  the  height 

of  the  barrier  and  a  is  the  hairier  half  width  at  the  zero  point  energy  position  in  the  reactant 
well.  This  expression  is  derived  from  the  WKB  approximation  to  the  wave  function  for 
particle  penetration  through  a  potential  barrier.  Assuming  that  this  barrier  is  modulated  by 
the  stretching  motion  between  the  naphthol  molecule  and  the  "ammonia  cluster"  simated 
near  the  -OH  group,  the  barrier  width  and  height  is  a  function  of  the  O— N  distance. 

Calculation  of  the  prottm  transfer  or  hairier  penetration  rate  as  a  function  of  the  van 
der  Waals  vibrational  excitation  in  the  cluster  requires  an  average  of  the  rate  constant  of  eq. 
(1)  over  the  wavefunction  for  the  "diatomic"  O  -N  distance.  Thus, 

k(v)  =  (Vvlk(a)fVv)  ,  (2) 

in  which  \\f^  is  taken  to  be  a  harmonic  oscillator  wavefunction. 

The  rate  constant  can  be  calculated  for  the  O  -N  vibrational  levels  of  cluster  Sj 

excitation  using  eq.  (2).  Note  that  the  vibrational  energy  in  the  cluster  is  never  sufficient  to 
excite  the  OH  stretch  either  directly  or  indirectly.  Calculation  of  requires  knowledge  of 

the  force  constant  and  reduced  mass  of  the  "diatomic"  O— N  stretch.  A  reduced  mass  of  38 
amu  is  calculated  for  the  O  -N  stretch  motion  with  naphthol  and  (NH^)^  taken  as  the  two 

moieties  in  the  "diatomic."  The  force  constant  is  taken  to  be  the  same  as  that  known  for  the 
phenol(NHj)|  system:^’*^  a  vibrational  energy  of  1 10  cm'*  is  calculated.  Several  other 

parameters  are  needed  to  evaluate  eq.  (2).  Specifically,  the  OH  stretch  in  Sj  is  chosen  as 
3000  cm'*  assuming  it  is  decreased  from  bare  phenol^**  by  hydrogen  bonding  and  the 
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0  --N  distance  in  Sj  is  taken  to  be  the  same  as  that  recently  determined  for  the  excited 
state  naphthol(NHj)j  complex,  2.60  The  barrier  halfwidth  at  the  zero  point  level  can 

be  estimated  for  the  classical  turning  points  of  the  harmonic  oscillator  to  be  ~  0.2A. 

This  model  has  only  one  adjustable  parameter;  the  barrier  height,  Ej^.  The 

vibrationless  proton  transfer  proton  transfer  rate  for  the  n  =  3  cluster  is  fit  to  the  observed 
rate  of  (57  ps)‘^  to  find  a  value  for  Ej^  of  65(X)  cm*^  The  transfer  rates  as  a  function  of  the 

vibrational  state  can  then  be  obtained.  The  details  of  the  above  calculations  are  given  in  the 
Appendix. 

Cluster  vibrational  excitation  in  Sj  is  generated  by  accessing  vibrational  states  of 


1-naphthol  which  are  then  assumed  to  undergo  rapid  relaxation  (intracluster  vibrational 

energy  redistribution  -  IVR)  to  the  cluster  van  der  Waals  modes.  Rapid  naphthol  to  van  der 
Waals  mode  IVR  for  the  1 -naphthol (Nlij)^^  clusters  is  a  good  assumption  based  on  the 

following  two  considerations:  high  cluster  density  of  vibrational  states;  and  known  IVR 
cluster  behavior  for  other  cluster  systems  (e.g.,  aniline(Ar)j,  (N2)j,  (CH^)j  and  indole 
(Ar)j^  and  (CH^jj  clusters).  The  energy  distribution  in  the  van  der  Waals  modes  is  then 
taken  to  be  statistical  and  can  be  readily  calculated.  The  probability  that  a  .specific  van  der 
Waals  mode  (in  this  case,  the  O  -N  "diatomic"  stretch)  is  populated  with  v  quanta  is  then 
calculated  using  the  Marcus  Rice  approximation^ 


in  which  E  is  the  vibrational  energy  in  the  cluster,  E^  is  the  energy  of  the  specific  van  der 


Waals  mode  and  S  is  the  number  of  van  der  Waals  modes  for  the  cluster.  The  products  of 
k(v)  and  for  the  modes  of  interest  are  summed  to  find  the  total  transfer  rate  for  a  given 
Sj  vibrational  energy  and 


k(E)  =  ^P^  k(v). 

V 


(4) 


Diatomic  harmor  ic  oscillator  levels  (hv  ~  1 10  cm'*)  to  v  =  3  are  found  to  contribute  to  the 
observed  rates  with  up  to  1500  cm'*  of  vibrational  energy  in  the  n  =  3  cluster.  Table  in 

presents  the  results  of  these  calculations:  agreement  between  the  calculated  transfer  rates 
(Xj'*)  based  on  the  above  model  and  the  experimental  results  is  excellent 

Isotopic  substitution  (H/D)  affects  the  proton  transfer  rate  in  Sj  for  the 
l  -naphthol(NH2)3  cluster  in  several  ways:  1.  the  OD  stretch  has  a  zero  point  energy  which 

is  much  smaller  than  that  of  the  OH  stretch  and  therefore  both  the  barrier  height  and  width 
increase;  2.  the  van  der  Waals  vibrational  reduced  mass  increases  for  (NDj)^  and  thus  the 

van  der  W  Js  vibrational  energy  decreases;  and  3.  the  proton  mass  is  doubled.  The  proton 
mass  appears  directly  in  the  exponential  (eq.  (1))  and  as  m  increases  the  rate  decreases. 
Increasing  the  barrier  height  Ej^  and  width  a  also  decreases  the  rate.  Decreasing  the  van  der 

Waals  vibrational  energy  also  decreases  the  transfer  rate  but  additionally  increases  the 
population  of  the  particular  mode.  Thus  the  deuteron  transfer  rate  for  the  cluster  should  be 
small  at  the  origin  and  should  increase  more  quickly  with  vibrational  energy  than  does  the 
proton  transfer  rate  for  the  protonated  cluster.  These  results  are  clearly  shown  in  Table  IE 
and  represent  excellent  qualitative  and  quantitative  agreement  between  the  model  and  the 

experimental  data.  Both  the  trends  and  the  details  are  well  described:  the  effect  of 
vibrational  energy  on  the  proton  transfer  rate  of  l-naphthoKNHj)^  clusters  is  reproduced 

and  the  decrease  of  the  isotope  effect  for  the  deuteron  transfer  rate  of  l-naphthol-djCNDj)^ 

clusters  as  a  function  of  vibraticMial  energy  is  predicted  as  well. 

Other  methods  of  calculating  the  O  -  N  separation  dependent  proton  transfer  rates, 
such  as  a  Fermi  Golden  Rule  approach,*^  can  also  be  employed.  In  this  approach,  two 
harmonic  potential  wells  (e.g.,  O-H  -N  and,  O- -H-N)  are  considered  to  be  coupled 
by  an  intermolecular  term  in  the  Hamiltonian.  Inclusion  of  the  van  der  Waals  modes  into 
this  approximation  involves  integration  of  the  coupling  term  over  the  proton  and  van  der 
Waals  mode  wavefimetions  for  all  initial  and  final  states  peculated  at  a  given  temperature  of 
the  system.  Such  a  procedure  requires  the  reaction  exothermicity  and  a  functional  fram  for 


the  variation  of  the  coupling  as  a  function  of  well  separation.  In  the  present  study  we 
employ  the  barrier  penetration  approach:  this  approach  is  calculationally  straightforward 
and  leads  to  a  clear  qualitative  physical  picture  of  the  proton  transfer  process. 

V.  CONCLUSIONS 

Excited  state  proton  rates  are  measured  and  modeled  for  l-naphthol(NHj)^,  n  =  2, 
3,  4,  5,  clusters  as  a  function  of  vibrational  energy  in  the  Sj  state  and  as  a  function  of 

hydrogen/deuterium  isotope  substitution.  These  rates  are  summarized  in  Tables  I,  II  and 
III. 

Proton  transfer  occurs  for  n  =  3  and  4  clusters  in  the  Sj  state  and  is  suggested  to 
occur  in  the  Sq  state  for  n  >  5.  No  transfer  is  found  for  n  <  3  clusters.  Substitution  of 
deuterium  for  hydrogen  in  these  clusters  (i.e.,  l-naphthol-djCND^)^)  has  dramatic  effects 
on  the  observed  transfer  rates.  At  the  origin  of  the  Sj4-  Sq  transition  the  kinetic  isotope 
effect  is  at  least  a  factor  of  20  (D^  transfer  is  slower  than  If'’  transfer)  but  at  an  energy  of 
1400  cm'^  the  kinetic  isotope  effect  is  about  a  factor  of  six.  The  proposed  model  for  proton 
transfer  in  these  clusters  equates  proton  transfer  with  a  simple  barrier  penetration  problem. 
The  height  and  width  of  the  barrier  are  modulated  by  vibrational  energy  in  the  cluster 

intermolecular  van  der  Waals  modes.  The  deuterium  isotope  effect  is  related  to  the  OD  zero 
point  energy,  the  reduced  van  der  Waals  mode  energy  for  the  l-naphthol-djCNDj)^  cluster, 

and  the  H/D  masses.  The  general  trends  and  detailed  rates  observed  arc  well  described  by 
the  model. 
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APPENDIX 


In  this  appendix  we  obtain  expressions  for  the  proton  transfer  barrier  height  and 
width  as  a  function  of  the  van  der  Waals  vibrational  displacements  and  derive  equations  for 
the  proton  barrier  crossing  rate  constant  as  a  function  of  the  van  der  Waals  vibrational  level. 

Assume  initially  that  the  surface  presented  in  Figure  3  consists  of  two  harmonic 
wells  (parabolas).  The  "barrier"  between  the  two  wells  is  represented  by  the  area  below  the 
intersection  of  the  two  parabolas,  but  above  the  zero  point  energy  of  the  O-H  vibration  in 
the  left-hand  well.  If  the  parabolas  are  symmetric  (i.e.,  the  vibrational  frequency  and  the 
minimum  energy  in  each  well  are  equal),  the  expression  for  the  barrier  height  at  the 
intersection  of  the  two  parabolas  as  a  function  of  the  separation  of  their  minima  is 

'  X  ^2 

E,  =  E„  (AI) 

V 

in  which  Eq  is  the  barrier  height  at  equilibrium  O  -N  separation  and  x  is  the  deviation  of 
this  separation  from  its  equilibrium  value.  The  variation  of  the  barrier  halfwidth  with  this 
deviation  is 


(A2) 


The  effect  of  small  asymmetries  of  the  size  expected  for  this  system  (^000  cm'^  of 
minimum  energy  difference  between  the  two  wells  and  frequency  differences  of  10%  to 
20%)  on  the  accuracy  of  eq.  (Al)  is  sufficiently  small  (<5%)  to  ignore.  Substituting  these 
intoeq.  (1)  gives 


k(x)  =  vexp 


(A3) 


Eq.  (A3)  may  now  be  substituted  into  eq.  (2)  and,  assuming  harmonic  oscillator 
wavefunctions  for  the  O  -  N  vibration,  the  integrals  which  must  then  be  performed  are 


of  the  form 
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v  =  1 
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v  =  3 


k  (V  =  3)  = 


yl 

3 


^^(0^/2 


exp 


(All) 
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Table  I.  Measured  Time  Constants  for  the  Nwi-Deuterated  Clusters  as  a  Function  of 
Vibrational  Energy. 


Mass 

n 

Excitation 

Energy  (cm-*) 

'Cl 

(ps) 

'C2 

(ps) 

3 

31,100  (0^ 

57(15) 

440(125) 

32,040 

43(11) 

500  (74) 

32,490 

12  (3) 

119(10) 

4 

30,800  (0^ 

71  (31) 

800  (400) 

32,040 

70  (30) 

600(100) 

32,490 

33  (20) 

1000  (500) 

Table  H.  Effect  of  Deuteration  on  Measured  Time  Constant  as  a  Function  of  Vibrational 
Energy. 


Mass 

H/D 

Excitation 

Energy  (cm-l) 

Xl 

(ps) 

X2 

(ps) 

H 

31,100 

57 

440 

D 

31,100 

>1000 

— 

H 

32,490 

12 

119 

D 

32,490 

75 

>1000 

Table  III.  Comparison  of  Barrier  Penetration  Model  Calculations  with  Measured  Proton 
Transfer  Time  Constants.  See  text  for  parameters  used. 


Mass  n, 

H/D 

Vibrational 
Energy  (cm"') 

calculated^ 

(ps) 

Xj  observed 
(ps) 

3.  H 

0 

57  fit 

57 

3,D 

0 

2600 

>1000 

3,  H 

940 

28 

43 

3,  H 

1390 

13 

12 

3,D 

1390 

200 

75 

\  =  6500  cm'‘,  a  =  0.2A,  E  (OH)  =  3000  cm’*.  E  (vdW)  =  1 10  cm'^ 


FIGURE  CAPTIONS 


Figure  1 
Figure  2 


Figure  3 


The  picosecond  laser  apparatus,  SPDL  =  synchronously-pumped  dye  laser. 
Time-resolved  1  +  1  ionization  signal  scan  of  the  l-naphthoKNHj)^ 

cluster.  Excitation  energy,  31,100  cm'*;  ionization  energy,  29,000  cm'*, 
o  =  data,  solid  line  =  fit.  The  fit  parameters  are  instrument  limited  rise 
followed  by  a  biexponential  decay  =  60  ps  and  X2  =  500  ps. 

Schematic  diagram  of  the  potential  model  used  for  proton  transfer  rate 
calculations.  Dashed  line  is  the  parabolic  barrier  using  parameters  fit  to 
origin  n  =  3  data.  Solid  curve  is  the  barrier  model  including  a  harmonic 

potential  well  with  an  OH  vibrational  energy  of  3000  cm"'.  The  well  is 
centered  at  1 .0  A.  The  other  parameters  for  this  model  are  Ej^  =  6500 

cm'*,  a  =  0.2  A,  and  the  van  der  Waals  stretch  energy  is  110  cm'*  for 


l-naphthoKNHj)^. 
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